Functional renormalization group approach to interacting bosons at zero temperature 
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We investigate the single-particle spectral density of interacting bosons within the non-perturbative 
functional renormalization group technique. The flow equations for a Bose gas are derived in a 
scheme which treats the two-particle density-density correlations exactly but neglects irreducible 
correlations among three and more particles. These flow equations are solved within a truncation 
which allows to extract the complete frequency and momentum structure of the normal and anoma- 
lous self-energies. Both the asymptotic small momentum regime, where perturbation regime fails, 
as well as the perturbative regime at larger momenta are well described within a single unified 
approach. The self-energies do not exhibit any infrared divergences, satisfy the U(l) symmetry con- 
straints, and are in accordance with the Nepomnyashchy relation which states that the anomalous 
self-energy vanishes at zero momentum and zero frequency. From the self-energies we extract the 
single-particle spectral density of the two-dimensional Bose gas. The dispersion is found to be of the 
Bogoliubov form and shows the crossover from linear Goldstone modes to the quadratic behavior of 
quasi-free bosons. The damping of the quasiparticles is found to be in accordance with the standard 
Beliaev damping. We furthermore recover the exact asymptotic limit of the propagators derived by 
Gavoret and Nozieres and discuss the nature of the non-analyticities of the self-energies in the very 
small momentum regime. 

PACS numbers: 05.10.Cc, 05.30. Jp, 03.75.Hh 



I. INTRODUCTION 

Interacting Bose gases in the continuum have been 
studied with various approaches for more than half a 
century, yet a controlled approach exists only for the 
weakly interacting variety and even there the asymp- 
totic regime at small frequencies and momenta proved 
to be non-perturbative. The first qualitatively correct 
description of the excitation spectrum of weakly interact- 
ing bosons was given by Bogoliubo\*i using a mean-field 
approximation. The Bogoliubov excitation spectrum ex- 
hibits the gapless linear Goldstone mode character at 
small wave vectors and approaches the quadratic disper- 
sion of free bosons as wave vectors become large and in- 
teraction effects become negligible. The crossover scale 
separating these two regimes is given by k c = 2mc where 
77i is the mass of the bosons and c the velocity of the 
Goldstone modes. In dimension D — 3, this picture has 
recently been confirmed experimentally using the Bragg 
spectroscopy technique on cold atoms £ 

Beliaev-2. went beyond the Bogoliubov approximation 
and calculated the self-energies of the quasiparticles to 
second order in the effective interaction, allowing both 
to extract corrections to the quasiparticle dispersion as 
well as to calculate their lifetime. However, as discussed 
in detail in Ref. perturbative approaches suffer from 
infrared (IR) divergences for all D < 3 which are diffi- 
cult to control at second order (in most physical quanti- 
ties the divergences cancel out) and essentially impossi- 
ble to control at any higher order of perturbation theory. 
These divergences have their origin in the U(l) symme- 
try of the model and can be traced to the divergence of 



the longitudinal correlation function^ - — Nepomnyashchy 
and Nepomnyashchy 5 showed rigorously that the IR di- 
vergences, together with Ward identities originating from 
the 17(1) symmetry of the model, lead to the surprising 
result that the anomalous self-energy £ (k,u)) vanishes 
at zero momentum and frequency. This result cannot 
be reproduced by finite order perturbation theory, and is 
in sharp contrast with both the Bogoliubov and the Be- 
liaev approach where the anomalous self-energy at zero 
momentum and zero frequency is finite. Furthermore, 
in the perturbative approaches a finite value of T, A (0) is 
a key quantity which controls both the crossover to the 
Goldstone regime as well as the velocity of the Goldstone 
modes. Nepomnyashchy and Nepomnyashchy showed 
further that the correct structure of the self-energies 
contains non-analytic terms which are dominant in the 
asymptotic limit of small frequencies and momenta. This 
non-perturbative structure arises for modes with mo- 
menta smaller than the generalized Ginzburg scale kc, 
which, for weakly interacting bosons, is much smaller 
than k c = 2mc. A correct analysis of these non-analytic 
terms allows to recover the Goldstone modes within a 
framework where X (0) = 0. There are a number of 
approaches which avoid the IR divergences, the earliest 
is due to Popov who separated the low energy modes 
from the high energy modes and treated the low energy 
sector within a hydrodynamic theory for the phase de- 
grees of freedom, which is free of IR divergences.— An- 
other approach is to artificially break the U(X) symme- 
try by a small symmetry breaking term which can in 
the end be removed to yield finite results for all physi- 
cal observables.— A natural framework to deal with IR 
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divergences is the renormalization group (RG). A field 
theoretical RG analysis of the interacting Bose gas at 
asymptotically small frequencies and momenta has been 
developed in Ref. @, which recovers the asymptotic be- 
havior of the self-energies and the single particle Green's 
functions. This required however a very careful analysis 
of Ward identities. 

In several recent publications the non-perturbative or 
functional renormalization group (FRG) has been ap- 
plied to the interacting Bose gas.— ~— The FRG is based 
on an exact flow equation of the generating functional of 
irreducible vertex functions (average effective action)^ 
Several approaches were based on a derivative expansion 
of the effective action which is sufficient to find a rela- 
tively simple description of the asymptotic regimei 12 i 13 
One advantage of the FRG approach is that the Ward 
identities associated with the U(l) symmetry are auto- 
matically obeyed if the truncation does not violate the 
invariance of the effective action. Another advantage of 
the FRG is that it can be employed to study the full mo- 
mentum and frequency dependence of the self-energies. 
See Refs. [l9l - l2l| for FRG calculations of the momentum 
dependent self-energy of classical models. In Ref. [l4| 
we showed how the FRG can be used to calculate the 
spectral density of the interacting Bose gas. This was 
based on a truncation of the exact FRG flow equation 
which respects the Nepomnyashchy relation and recov- 
ers the non-analytic structure of the low energy sector. 
Moreover, our FRG approach also describes the full dis- 
persion of the Bogoliubov spectrum as well as the Beliaev 
damping of the quasiparticles. This is the first approach 
which succeeds to describe both the non-perturbative as- 
pects which appear at the very small wavevector scale 
k <^kc, as well as the perturbative ones, which are con- 
trolled by the larger scale k c = 2mc, within a unified 
framework. Recently Dupui o 15 ! 16 developed a similar ap- 
proach and clarified the role of ka in the flow equations, 
but did not resolve the Beliaev damping of the quasipar- 
ticles. 

Here we present in detail our approach of Ref. [3| and 
investigate analytically the different scaling regimes of 
the flow equations. We further derive the flow equa- 
tions for interacting Bose systems which treat the ir- 
reducible density-density interactions exactly and which 
has not been published previously. This set of equations, 
while difficult to solve, would potentially allow to ana- 
lyze strongly interacting yet dilute Bose systems, which 
can now be realized in experiments and whose spectrum 
shows clear deviations from the prediction of the Beliaev 
approach. 22 

This work consists of the following parts: In Sec. [II] we 
define the model and briefly review the structure of the 
one-particle Green's functions, i.e. discuss the excitations 
and their damping. We further summarize the structure 
of the self-energies in the non-perturbative regime. In 
Sec. Mil we discuss the FRG approach, introduce our ap- 
proximation and discuss both the general structure of the 
flow equations as well as the results of a derivative ex- 



pansion in subsection IIII Al for which we present results 
in D = 2. In Sec. IIVI we finally discuss the results of a 
non-self-consistent solution of the flow equations which 
yields the full momentum and frequency dependence of 
the self-energies and allows to extract the single-particle 
spectral function. Numerical results are shown for D = 2. 
In Appendix [X] we briefly review the arguments from 
Ref. Q which lead to £^(0) = 0. An analysis of the 
derivative expansion is presented in Appendix [B] where 
the flow of the effective interaction is analyzed analyti- 
cally and the standard T-matrix results from diagram- 
matic approaches are recovered for general dimension D. 
Furthermore, in Appendix [C] we show how the Ginzburg 
scale fee emerges from the flow equations. In Appendix [Dl 
we further show analytically, how the exact results for the 
asymptotic form of the normal and anomalous propaga- 
tors is obtained from the non-self-consistent approach of 
Sec. IIVI We conclude our work with a summary in Sec.fVl 

II. SELF-ENERGIES OF THE INTERACTING 
BOSE GAS 

We consider a system of bosons at zero temperature 
with a repulsive interaction potential u A (fc). The bare 
microscopic action is defined in terms of the following 
functional of a complex field ip, 

S[i>,ilA, = - J 4> K ( iuj - £ fc + m)V>k 

K 

+ \ J u A ( k )P-KPK> (2- 1 ) 
K 

where K = (fe, iw) and k is a D-dimensional momentum 
variable with the absolute value k and iu> a bosonic Mat- 
subara frequency. At zero temperature the integration 
over K is defined as 

K 

The first term of Eq. (|2.II) describes non-interacting 
bosons with mass m and the dispersion e& = k 2 /2m. 
The chemical potential is denoted by fi. The second part 
represents the interaction where p K denotes the Fourier 
transform of the local density 

PK = J i>QipQ+K ■ (2.3) 
Q 

The subscript A of u Aa (k) indicates that the theory is 
assumed to be regularized in the ultraviolet (UV) by a 
finite cutoff Ao . This cutoff is related to the finite range 
of the interaction; for example, for hard core interactions 
Ao ~ 2ir/d, where d is the diameter of the particles. The 
U(l) symmetry of the model Eq. (|2.ip is spontaneously 
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broken in the ground state and the field ip acquires a 
finite grand canonical expectation value, 



0° = (^(r,r))^O. 



(2.4a) 



A macroscopic number of particles in the ground state 
forms a condensate with density p° = |0°| 2 . Because 
(f>° ^ 0, the Green's function has both normal and anoma- 
lous components. The Dyson equation for the Green's 
functions in the symmetry broken state can be compactly 
written in a 2 x 2 matrix in the space spanned by the field 
types ipK and as 



G- 1 {K) = G-\K)-V{K), 



with 



G{K) 



G N (K) 
G A (Kf 

Y, A {K)* 





G A {K) 
G N (-K) 

T, N (-K) 




and 



Go 1 (K)=iu-e k - 



This leads to 



G(K) = 

1 f -Go 1 (-K) + E If (-K) 
T, A (K)* 



V{K) 

where the denominator T>(K) is given by 



T, A {K) 
Go\K)^^ N 



(2.5) 

(2.6) 
(2.7) 
(2.8) 

(2.9) 
(2.10) 

^ {K) 



V{K) = -[G^(K)-E N (K)][G^(-K)-E 
+ \Z A (K)\ 2 . 



N , 



-K)] 
(2.11) 



We shall refer to G N (K) as the normal propagator and 
G A (K) as the anomalous propagator^ The U(l) sym- 
metry imposes constraints on the vertices among which 
the Hugenholtz-Pines relation for the self-energies 2 ^ 



£ iv (0)-5T(0) = /i 



(2.12) 



is the best known. This equation is obeyed by the Bogoli- 
ubov approach 1 which approximates both self-energies by 
the frequency independent leading order result E w {K) ~ 
P°[ua o (0) +UA (k)} and Z A o (K) sa p Q u Ao (k), while the 
condensate density is given by p° ~ m/ u a (0) ■ This yields 
undamped excitations with dispersion 



E, 



(2.13) 



For small k the dispersion has a Goldstone mode charac- 
ter with Ek ~ cq A; with the velocity 



Co 



(2.14) 



while at large momenta E k approaches the free particle 
dispersion e k — k 2 /2m. If UA (k) — u\ is independent 
of k, we can write 



E k 



where 



c fc^/l + k 2 /k t 



ko = 2mco 



(2.15) 



(2.16) 



is the mean-field estimate for the characteristic scale k c — 
2mc of the crossover from Goldstone modes to quasi-free 
bosons. 

The leading (second order in the interaction) many- 
body corrections to the Bogoliubov approximation have 
been calculated by Beliaev.— At this level of approxi- 
mation the single-particle excitations are damped. For 
k — > the damping j k has in three dimensions the forroS.^ 



(£>=3) 

Tk 



3fc 5 



6407rmp° 



while in two dimensions 25 - 26 



{D=2) 

% 



V3c Q k 3 



(2.17) 



(2.18) 



32irp° 

The general result in D dimensions can be written as^ 



where 



P + l 

a = 3 ' Kd- 



(d) 
Tk 



, 3-D 



32mp° 



„2D-1 



dx : 



-\l-x) 



D-l 



with 



2 d - 1 tt d / 2 T[D/2} 



(2.19) 



(2.20) 



(2.21) 



However, as discussed in detail in Ref. [J], already at sec- 
ond order the diagrams for the self-energy are in fact IR 
divergent and the divergent terms have to be separated 
from the non-divergent terms in order to obtain physi- 
cal results. These divergences are in fact related to an- 
other problem common to all perturbative approaches, 
which is that they violate an exact result obtained by 
Nepomnyashchy and Nepomnyashchy^ 5 who showed that 
the momentum and frequency independent part of the 
anomalous self-energy vanishes for D < 3, 



T, A (K = 0) = 0. 



(2.22) 



This result is a direct consequence of the IR divergences 
and a Ward identity relating the three point vertices 
to the two point vertices, see Appendix [S] While in 
the Bogoliubov and Beliaev approach a finite value of 
Ti A (K = 0) ensures the Goldstone character of the long 
wavelength modes, in the exact theory the Goldstone 
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modes are recovered by a non-analytic structure of the 
self-energies which have at small frequencies and mo- 
menta the form 5 - 

T, N (K) fa p°u{K) + p + iuj + a N e k + b N cu 2 , (2.23a) 



E A (AT) ps p°u{K) + a A e k + b A Lo 2 



(2.23b) 



where u(K) is a non-analytic function of u and fc and 
the coefficients ajv, &jv, a, a, b A depend on the interaction 
(note that the coefficient of the iu term is exactly equal 
to unity 5 ). Thus, while each of the self-energies is non- 
analytic for D < 3, their difference is analytic, so that 
the quantity 



a(K) = Y, N (K)-Y 1 A {K)~ p 



(2.24) 



has the expansion 

a(K) ss iuo + (ojv — a A )ek + {bN — b A )uj 2 ■ (2.25) 

We will use in our FRG approach below a scheme 
where this structure of the self-energies appears nat- 
urally. As long as all irreducible correlations beyond 
density-density correlations can be neglected, u{K) is 
in fact the fully renormalized density-density interaction 
which also completely determines the anomalous self- 
energy, i.e. Y, A (K) = p°u(K). 



III. FUNCTIONAL RENORMALIZATION 
GROUP APPROACH TO THE INTERACTING 
BOSE GAS 

We shall now set up the FRG equations for the inter- 
acting Bose gas . For recent reviews of the FRG method 
see Refs. [27H33|. The basic idea is simple and consists of 
introducing an IR cutoff A to regulate the theory which 
is then sent to zero in infmitesimally small steps. The 
FRG follows the evolution of the various vertex func- 
tions as the IR cutoff A is lowered and yields the true 
vertex functions for A — > 0. This will allow us to extract 
the excitation spectrum and the damping of quasiparti- 
cles from the renormalized self-energies within a theory 
which is free of IR divergences. 

The IR divergences which plague the perturbative 
approaches^ are removed if we regulate the theory by 
introducing a momentum cutoff A in the free propagator, 



Gq \(K) =iu- e k + p- R A (k) . 



(3.1) 



Here, R\(k) is a regulator function which removes the 
IR divergences arising from modes with k < A and will 
be specified later. The FRG approach is based on the 
cutoff-dependent effective action which is defined by 

r A [0, 4>] = C A [4>, <j>] - [ K GvX(K)<t>K , (3.2) 



A" 



where Ca[4>,4>] is the cutoff-dependent Legendre 
transform of the generating functional of connected 



Green's functions. The effective action is the 
generating functional of irreducible vertex functions, 
i.e. it generates cutoff-dependent irreducible vertices 



, K\) which in the condensed 



phase are defined via the functional Taylor expansion of 
the corresponding generating functional in powers of the 
fluctuations 84>k — ~ $K 0*1 



E 



i 



n\m\ 



K' 



KL K„ 



n,m— 

x ^A;+...+A;,A m +...+A! 
x T^ m \K[,...,K' n ,K n 
x S(j) K i . . . S(f>K'Scj) K . . . 



Ki 

h Kl ■ (3.3) 



We normalize T&[<fi,(p] such that, at the initial RG scale 
Ao , it reduces to the bare interaction minus the chemical 
potential term \i J K (f) K cj) K ^ Note that 0° now also has 
an explicit A-dependence which will be chosen such that 
r^ 1,0 - 1 and T 1 ^' 1 ^ vanish for all A. 34 For our purpose it 
is sufficient to approximate the effective action by the 
following ansatz, 



A 



A 



where 



Sp 



A 



(3.4) 



(3.5) 



is the Fourier-transform of the condensate density fluc- 
tuation p(X) — p A = 4>x4>x — Pa- Here, in analogy with 
the definition of K, we use X — (r. r), where r is a D- 
dimensional real space coordinate and r is the imaginary 
time. At the initial UV cutoff scale A = A Eq. I|3.4p 
is exact and ua (A) is just the bare interaction which 
enters Eq. (12. ip whereas a a vanishes initially, 



<7A W=0. 



(3.6) 



While the effective action (|3.4j) describes arbitrarily 
strong density-density interactions, processes which in- 
volve three or more particles are neglected. Thus, the 
ansatz is not expected to be accurate for dense systems 
but it can describe dilute yet strongly interacting ones. 
Note that the effective action is completely parameterized 
in terms of the scalar p A and the two cutoff-dependent 
scalar functions <ta(A) and u\(K). Eq. (|3.4[) represents 
a non-local potential approximation of the action and is 
explicitly U(l) invariant for any value of A even if the 
condensate density p A is finite. In the usual field ex- 
pansion, the U(l) symmetry of the model leads to Ward 
identities which relate higher order irreducible vertices 
to lower ones. In our approach the irreducible vertices 
are derived from an explicitly invariant effective action 
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which guarantees that all Ward identities are automati- 
cally obeyed. The normal self-energy 



and the anomalous self-energy 

Z£(K) = r£°\-K, K) = Yf 2) (K, -K) 



(3.7) 



(3.8) 



have the form (here and below we shall assume, without 
loss of generality, a real valued condensate wave function 



b A for notational convenience) 



(K) 



fj, + a A (K) 
P>a(K) . 



P>a(K), 



(3.9a) 
(3.9b) 



Note that Eqs. (|3.9al3.9b[) have the same structure as 
Eqs. (|2.23al2.23bp and we can already anticipate that 
ua(K) will become non- analytic for A — > while 
cr\(K) will remain analytic. Moreover, since cr A (0) van- 
ishes, Eqs. (|3.9al3.9bp obey the Hugenholtz-Pines rela- 
tion (l2"T2")l for all A. Note that Eq. has no terms 
linear in the field 8<f> which for A = Ao is achieved by 
fixing the initial condensate density as 



where 



u A n = u A (0). 



(3.10) 



(3.11) 



Eq. (|3.4p also fixes the three- and four-point vertices of 
the theory which have the symmetrized form 

T^iK'^K^K,) = <[u A (^) + <x A (i^)], (3.12) 



r^ 2) (^;K 2 ,K!) = cbHuAiK^ + u A (K 2 )}, (3.13) 

(3.14) 

These vertices are not independent but are fixed entirely 
by <p A and u A (K) which can be extracted from the self- 
energies. This dependence is again a consequence of the 
U(l) symmetry of the model. 

For the FRG we also need the cutoff-dependent full 
propagators which can again be written in the form of a 
2 x 2— matrix in field type spaced 



where 



Ga{K) 

G N a{K) 
Gi(K) 



G N K {K) Gf(K) 
g a( k T Ga(~K) 



-G- K U-K) + ^{-K) 



Va{K) 



■Da(K) 



(3.15) 

(3.16a) 
(3.16b) 



and the denominator T> A (K) is given by the expression 

V A (K) = [w + e k -n + HZ(-K)] 

x [iw - e fc + a* - Sa (K)} + \^i(K)\ 2 . (3.17) 

The exact flow equations for the effective action T A [4>, <fi] 
has the forrn^ (here we do not keep track of the flow of 
the constant part of T A [<f>, <fr], i.e. the free energy^) 



d A r A [<M = ^Tv{d A G-\[T^ - G-XV 1 } 



(3.18) 



where 



,(2)r 



/ S 2 T 
\ 



5 2 T\ 



6 2 T\(j>,, 



5 2 T\ 



\ 
I 



(3.19) 



Eq. (|3.18j) . in conjunction with the field expansion (|3.3|) . 
yields flow equations for the irreducible vertex functions, 
see e.g. Refs. [Hll for details. 

To ensure that the field expansion p. 31) is always an 
expansion around the true minimum (j> A , we require that 
r^ 1,0 -* and r^' 1 - 1 vanish for all A. Since for A = Ao this 
is already achieved by the mean field Eq. (|3.10[) . we will 

only need to enforce c^r^' 1 ^ = SaF^ 1,0 ' = 0. This yields 
the flow of the order parameter p A , 



dA P °A = 



1 



ua(0) 



Gf(Q)[u A (Q)+u A (0)] 
+Gi(Q)u A (Q)]} , (3.20) 



which is shown diagrammatically in Fig. [TJ The single 
scale propagators G A and G A which appear in Eq. (|3.20[) 
are defined via the matrix equation 



G<lh) ] - -G A (K)[8 A GoX{K)]G A {K) . 

(3.21) 



Gi(K)* Gf(-K) 
where 



G X(K) 
G Q X(-K) 



(3.22) 



The exact FRG flow equations for the normal and anoma- 
lous self-energies are shown diagrammatically in Figs. [5] 
and [3] Using the relations Eqs. (13. 12113. 14|) for the three- 
and four-point vertices, the flow equations for the self- 
energies reduce to 
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-Pa y G^(Q){g^(K - Q)[u A (Q) + u A (K - Q)] 2 + G^(Q - K)[u A (K) + u A (K - Q)] 2 
Q 

+Gl (Q + K) [u A {Q) + u A {K)] 2 + 2Gi(K - Q) [u A (Q) +u A {K- Q)] [u A (K) + u A (K - Q)] } 



-2pl J Gi(Q){Gi(K-Q)[u A (K) + u A (K-Q)]+G%(K-Q)[u A (Q)+u A (K-Q)]} 
Q 

x [u A {K) + u A (Q)] , (3.23a) 



d A x£(K) 



-Pa 



u A {K) 
ua(0) 



G 



N 



(Q) K(Q) + u A (o)] + g£(Q)u a (q)} -I I G£(Q) [u a (K + Q) + u a (K - Q)] 



a y G% (Q){ [G A (Q-K) + G% (K - Q)) [u A {K) + u A (K - Q)] [u A {K) + u A (Q)] 
Q 

+ (G£(K + Q) [u A (Q) + u A (K + Q)] + G A (K + Q) [u A {Q) + u A (K + Q)] ) [u A {Q) + u A (K)} } 

Pa J Gi(Q){Gi(K - Q)([u A (Q) + u A (K - Q)} 2 + [u A (K) + u A (K - Q)] 2 ) + Gi{K + Q) [u A {K) + u A (Q)] 2 
Q 

+ [G A (K -Q) + G a (Q - K)\ [u A {K) + u A (K- Q)) [u A {Q) +u A (K~ Q)] } , (3.23b) 



which, together with Eqs. (|3.9al3.9bj) . yield the flow 
equations for a A (K) and u A {K). We stress that up 
to this point the only approximation we have made is 
that we only kept those irreducible correlations which 
are explicitly stated in the Eq. (13.41) for Ta[$, <f>], i.e. 
we only consider two-body density-density correlations. 
However, these can be arbitrarily strong and the set 
of Eqs. q3.2QI3.23al3.23bj) allow to calculate p° A , u A (K), 
and a A (K), which completely determine all vertices of 
Eq. (I3.4[) . without any further approximation. In this 
respect, our approach goes be yon d the set of equations 
given in Appendix C of Ref. [16( where additional ap- 
proximations were made. While of similar complexity, 
the equations derived by Dupuisi^ do not treat the full 
frequency- and momentum-dependence of the three- and 
four-point vertices, which appear in the flow of the self- 
energies, exactly. In some diagrams the dependence of 
these vertices on the internal frequency and momentum 
is neglected. 



Eqs. (|3.20l3.23al3.23b[) form a set of integro-differential 
equations which in principle must be solved self- 
consistently. While similar flow equations were re- 
cently solved exactly for a model of classical crystalline 



FIG. 1: (Color online) Diagrammatic representation of the 
exact FRG flow equation for the condensate density. The 
circles with two external legs on the left-hand side denote 
the normal and anomalous self-energy. The circles with three 
external legs on the right-hand side denote irreducible three 
point vertices. Crossed circles with a dot denote the deriva- 
tive of the order parameter with respect to the RG cutoff A, 
and crossed solid arrows represent normal or anomalous single 
scale propagators as defined in Eq. (|3.21[) . Arrows pointing 
out of a vertex represent the fields tp, while incoming arrows 
represent fields tp. 



membranes,— for the present quantum system the cal- 
culation of an exact solution is an extremely difficult 
problem which we will not attempt to solve. Rather, 
we will use a non-self-consistent approximation which is 
motivated by the structure of perturbation theory. Note 
that second order perturbation theory (including its di- 
vergences) is recovered if we replace on the right hand 
side of Eqs. p.2QI3.23al3.23b|) the full function u A (K) 
by the constant u Ao and use a A = 0. 



Before we calculate the full momentum and frequency 
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FIG. 2: (Color online) Diagrammatic representation of the 
exact FRG flow equation for the normal self-energy S^v {K). 
The circle with four external legs denotes the irreducible four- 
point vertex, while arrows represent cutoff-dependent normal 
or anomalous propagators. The dot above the two-point ver- 
tex on the left-hand side represents a derivative with respect 
to A. The other symbols are defined as in Fig. [1] 




FIG. 3: (Color online) Diagrammatic representation of the ex- 
act FRG flow equation for the anomalous self-energy (K) . 
All symbols are defined as in Figs. [1] and [2] 



d A X»(K) 



We now further employ a derivative expansion of cr A (K) , 
which is analytic also for A — > 0, where we keep only the 
leading order terms in an expansion in u) and k. We thus 



dependence of the self-energies using a non-self-consistent 
scheme, we first consider a derivative expansion of T\ in 
the next subsection. 



A. Derivative expansion 



A derivative expansion of u A {K) is problematic, since 
ua(K) becomes non-analytic for A — > 0. We therefore 
approximate it only by its constant part, 



ua(-K")w«a(0) = «a. (3-24) 



This leads to a great simplification for the three- and 
four-point vertices. With the approximation (I3.24[) one 
finds from Eq. (|3.20p the simplified form for the flow of 
the condensate 17,21 



d A p° A = J (2G%(K) + G£(Kj), (3.25) 

K 



and from Eqs. (|3.23a|3.23b|) and (|3~25j) the simplified 
flows for the self-energies , 14 > 17 



approximate 

a A {K) a iw (1 - Y A ) + e k {Z^ - 1) + u 2 V A , (3.27) 



2u A j {G^(Q) + Gi(Q)}-4ulp° A j {G^(Q)[G^(Q + K) + G^(Q-K) + G^(-Q + K) 

Q Q 
+2G£(Q - K)] + 2Gi{Q) [Gf{Q + K) + G^(Q + K)] } , (3.26a) 

2u A J G^ (Q) ~ \u\p\ J [g£ (Q) [G% (Q + K) + G N h (Q-K) + G£(Q + K) + G A (Q - K)] 
Q Q 

+Gi{Q) [G A (Q + K) + G% (Q-K) + 3G^(Q + K)] } . (3.26b) 
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where initially we have Za — Y\ — 1 and Va = 0. 
With this truncation, the normal and anomalous propa- 
gators are simply 



G%{K) 



Gf(K) 



Y A ioj + Zl l e k + A A + R A (k) + V a oj 2 



(3.28a) 
(3.28b) 



where 

V A {K) = Y^ 2 + [V A u 2 + Z A 1 e k + R A {k)] 

x [2A A + V a lo 2 + Z^e k + R A (k)] , (3.29) 

and we defined 

A A = P ° A u A . (3.30) 

The corresponding single-scale propagators follow from 
Eq. (|3.2f [) . We now employ the Litim regulator— defined 
through 

R A (k) = (1 - 4.o) (2m2 A )- 1 (A 2 - k 2 ) 6 (A 2 - k 2 ) , 

(3.31) 

which leads to a simplified form of the flow equations 
since the integration over the internal momentum can be 
performed trivially, 

rI D k A D+1 
^d A R A {k)F{k 2 ) = KA J-(A 2 ), (3.32) 

for any function jF(k 2 ). Here, we defined 

K A = K D [l-r) A /(D + 2)]/D, (3.33) 

where 

jf A = A8 A lnZ A (3.34) 

is the scaling dimension of Z A and Kjj is defined in 
Eq. (j2~2i"T) . From Eqs. (|3.25l3.26bj) wc find the flow of 
p A and ma 



dAP°A = 4 



a d+i k A r du 



9^ 2^, 



3 Jp), ,2r, 



2mZ A J 2tt ^ D\(iuj) 



d A u A = Au\ 



A d+1 k a [ du) 



3 „(«), ,2r, 



2mZ A J 2^^D\{iu) 



(3.35a) 



-, (3.35b) 



where 



D A (ioj) = V A (A,iw) 

= Y 2 A u? + [e A + V A uj 2 ][e A + V A u 2 + 2A A ], (3.36) 

and we used the shorthand e A = e A /Z A . The coefficients 
entering the flow (|3.35ap are 



= A + 6aA a + A 2 : 

'1 

» 



Va(2ca + Aa) — Y A , 



(3.37a) 
(3.37b) 
(3.37c) 



and those entering (|3.35b[) are 



» 
-o 

» 



5e A + 3? A A A + 6? A A A + 4A A , (3.38a) 



cr = W A (5ei + 2e A A A + 2A|) 

-F A 2 (7A A + lle A ) , (3.38b) 

4 M) = V A [W A (5e A + A A ) - 11F A 2 ] , (3.38c) 

4 U) = 5V A 3 . 



(3.38d) 



To find the flow of the parameters entering cja{K), we 
must expand the flow of £ A (K) , (I3.26ap , to second order 
in uj and k. This yields 



9aYa = 



2rnZ A J 2w£^Dl(i 



d A Z A = 4p° A u 2 A 



A D+1 K n f du 1 



2mD J 2tt D 2 A (iu) 



d A V A = 8p A < 



2 A D+1 KA /"dw^c^w 2 " 



2mZ A 7 2ir f^ Q Dl{w) 



, (3.39a) 
(3.39b) 
(3.39c) 



where the coefficients entering the flow of 1a are 



Sv) - ?2 



= ei — 2e A A A — 2A 



2 

A • 



Y 2 + 2(e A -A A )V A , 



(3.40a) 
(3.40b) 

VI , (3.40c) 
and the coefficients entering the flow of V A are 



M - T/2 



» _ 

'1 



-y A 2 (e~ A + A A ) - e A (e A + 2A A )V A ,(3.41a) 



2V A [Y 2 + V A (i A + A A )} 
3K 3 . 



(3.41b) 
(3.41c) 



Before we analyze the results of the derivative expansion, 
let us briefly consider the scaling dimensions of the cou- 
pling parameters entering the theory. If we define the 
momentum to have scale one, [k] = 1, and the dimension 
of frequency to be equal to the cutoff dependent dynam- 
ical exponent [lo] — z A , we find 

[Va] = 2-2z A -77 A , (3.42a) 
[ua] = 4-D-z a -2 Va , (3.42b) 
D - 2 + z A + r] z A , (3.42c) 



[p°A 



where the scaling dimensions of Z A and V A are defined 
by 



[Za] = Va = Ad A Z A , 
[Y A ] = vl = A9 A y A , 

which fixes the dynamical exponent, 
z\ = 2 - Tj A - rf A . 



(3.43a) 
(3.43b) 



(3.44) 



Initially, we have ij Ao = rf Aa = and z Ao = 2 which is the 
correct dynamical exponent for non-interacting bosons. 
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It is then useful to give the bare interaction strength ua 
and the chemical potential in dimensionless form, 



D-1 



Uq = 2mw Ao A 
fl = 2m^A Ao 2 



(3.45a) 
(3.45b) 



Note that at T = 0, the scaling dimension rf K plays no 
important role since it vanishes at A = 0. However, at the 
finite temperature transition of the Bose gas the critical 
fluctuations lead to a finite rf = limA^o r /A- This was 
analyzed in Ref. [13] where the anomalous dimension rf 
was evaluated at the critical point. The relevant scaling 
dimension at T — is however 77^ which controls also the 
crossover from the za = 2 regime to the Z\ = 1 regime 
which is characteristic of the Goldstone modes . 12 ' 13 
Equations similar to Eqs. (j3"35al l3~35bl l33i?aTl339"c]) 



were already discussed in Refs. 
summarize the results for 2 < D < 3 
densate density p A flows to some finite limit p^-j-o = 
0, while the coupling constant ma vanishes for A 



HEf and we briefly 
At T = the con- 

P > 
■> 



which ensures 



(K = 0) = n and T,^ (K = 0) = 0, 
in accordance with both the Hugenholtz-Pines relation 2 ^ 
and the Nepomnyashchy identity. Furthermore, Y A van- 
ishes for A — ► 0, again in accordance with exact results. 5 



Since V A and Za are finite for A 



0, the low en- 



ergy modes are phonons with a linear dispersion and 
velocityi2r-i£ 



1 



V2mVZ 



(3.46) 



with V = Va-j.0 and Z = Zj\^o- In Fig. 0] and Fig. [5] we 
show for D — 2 the flows for the coupling parameters and 



for the dynamical exponent z. 



v 
Va 



rj A which il- 



lustrates the crossover from the free particle regime with 
2A = 2 to the Goldstone regime, where za = 1- As 
pointed out by Dupuis in Refs. flBTflrjl ]. the crossover of 
za is governed by the generalized Ginzburg scale kc, 
which is indicated in Fig. 2] and Fig. [5] and is extracted 
from the flow equations in Appendix [Cj Also indicated 
is the crossover scale k c = 2mc. Evidently, kc is the rele- 
vant scale controlling the flow of the dynamical exponent 
za ' 15]1G On the other hand, the scale which determines 
the crossover from a linear to a quadratic behavior in the 
full quasiparticle dispersion is not kc but the crossover 
scale k c , as is expected from standard perturbative ap- 
proaches and as we will show in the next section when 
we calculate the full frequency and momentum depen- 
dence of the self- energies. It is perhaps surprising that 
the crossover in the momentum dependence of the quasi- 
particle dispersion happens at a completely different scale 
as the crossover in the A-dependence of the parameters 
entering the derivative expansion. In the weak-coupling 
regime this is however easily understood, since a quasi- 
particle with momentum k in the range kc *C k <C k c 
feels the largely unrenormalized bare value ma of the 
interaction potential which is sufficient to render its dis- 
persion linear, whereas the flow of the parameters en- 
tering the derivative expansion are sensitive only to the 




J ln(A /k c ) 2 ln(A /k G ) 4 

ln(A /A) 



a. 



°< 0.95 




ln(A /k c ) 2 ln(A /k G ) 4 

ln(A /A) 



FIG. 4: Typical RG flows of the interaction parameter u A 
and the condensate density p A for D = 2. The initial val- 
ues are /to = 2m/iA^" 2 = 0.4 and Uo ~ 2muA A^ -2 = 4. 
The Ginzburg scale &g and the crossover scale k c — 2mc are 
indicated by vertical lines. 



non-perturbative effects associated with the emergence 
of non-analyticity at scales of the order of the Ginzburg 
scale kc- 

The lowest order derivative expansion is thus incapable 
of describing the quasiparticle dispersion beyond the 
asymptotic limit. Strictly speaking, even this limit is not 
correctly reproduced since the derivative expansion yields 
the unphysical result limA^o ^a(K) ~ fi m A->o /°a u a = 
for all K, and as a consequence limA^o Gj^{K) —> 0. This 
violates an exact result by Gavoret and Nozieres 3 - 7 which 
is given in Eq. (|A1|) . Nonetheless, it is possible to extract 
the correct asymptotic form of the propagators if one cuts 
off the flow at a finite A and uses A 
Refs. dHHi]. 



c 2 k 2 , see 



IV. EXCITATION ENERGY AND DAMPING 
OF QUASIPARTICLES 

We now turn to the calculation of the complete fre- 
quency and momentum dependence of the self-energies 
from the FRG equations. As pointed out in Sec. IIII1 a 
completely self-consistent solution of the flow equations 
(13.201 13.23al I3.23b|) is numerically extremely difficult and 
we therefore use an approximation in which we treat 
only the coupling parameters of the derivative expan- 
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1.2 



0.8 



0.4 







V A /V 




;-- z, 










ln(A /k c ) 2 ln(A /k G ) 4 


6 


i 



ln(A /A) 



2.0, 



1.5 



1.0 



0.5 



A 2-ll A -r\ 



ln(A /k c ) 2 ln(A /k G ) 4 6 8 

ln(A /A) 

FIG. 5: Typical RG flows of the coupling parameters Y A , 
Z A , V A and the dynamical exponent z A and both anomalous 
dimensions rf A and rj A in D — 2. The initial values are the 
same as in Fig. 2] 



sion of Sec. MI Al self-consistentlv. while all higher powers 
of frequency and momentum are determined in a non- 
self-consistent manner. A similar truncation has previ- 
ously been tested for the symmetry-broken phase of a </> 4 
theory.— A non-self-consistent solution can be quite sim- 
ply achieved, starting from the approximate flow equa- 
tions given in Eqs. fl3.26al I3.26b[) . On the right hand 
side of these equations we use the approximate form of 
the propagators given in Eqs. (I3.28al I3.28b|) which are 
based on the derivative expansion, see Eq. (I3.24j) and 
Eq. (|3.27[) . The same approximation is employed for the 
single scale propagators. Note that Eqs. Q3.28a[ I3.28b|) 
have a sensible large K behavior which for any finite 
A and all k > A arc similar to the standard first or- 
der result of the Bogoliubov approximation but with a 
renormalized interaction parameter u\, a renormalized 
condensate density p A , a mass renormalization expressed 
by Z\ and further a renormalized frequency dependence 
as expressed through Y\ and Va- 

Using again the Litim regulator (I3.31[) . the flows of 
the coupling parameters ma, Pa, Va, Z\ and Y A can be 
simply obtained from the solutions of Eqs. (|3.35al I3.35bl 
I3.39aj|3"!39cl) . With these approximations, all parameters 
entering the right hand side of Eqs. Q3.26a[ I3.26bj) arc 
completely determined. For any K, we can thus simply 
perform the integration over A. With appropriate bound- 
ary conditions this yields both T, N (K) and T, A (K). Be- 



cause of the structural similarity of both the diagrams 
and the propagators in this approximation to the Beli- 
aev theory, we expect to reproduce qualitatively Beliaev's 
perturbative results for both the spectrum and the damp- 
ing, bar their divergences. However, since already the 
derivative expansion contains the key information about 
the non-analytic structure of the self-energies, our ap- 
proach also captures these, as discussed in more detail in 
Appendix |DJ 

All information about the quasiparticle properties is 
encoded in the single-particle spectral density which is 
related to the imaginary part of the real-frequency nor- 
mal Green's function, 



A{k,u) = -2ImG iV {k,iw 



•tO). 



(4.1) 



Its calculation requires an analytical continuation to real 
frequencies. Here we apply the standard Pade approx- 
imant technique^ which has the advantage that it can 
be easily implemented numerically. Although it has been 
originally proposed for dealing with systems at finite tem- 
peratures, it can also be used for zero-temperature calcu- 
lation if the input functions are sufficiently well resolved. 
If it does converge, the technique is furthermore quite 
accurate. Applied to our problem it proved to be very 
robust and the results for the spectral density quickly 
converged for momenta which are not too small. For 
the actual data used to calculate the spectral density we 
used 450 Matsubara frequency data points for each mo- 
mentum k. Even though the damping is quite small, this 
allowed to accurately extract the width of the resonance 
and the damping. 

The first task is to calculate the Matsubara normal 
Green's function with satisfactory resolution. A typical 
result for the normal Matsubara Green's function cal- 
culated within the non-self-consistent FRG scheme just 
described is shown in Fig. [6j The initial conditions for 
the condensate and interaction were chosen in such a way 
that the ratio k c /k G remains sufficiently large (see leg- 
end below Fig. We further kept the ratio c/c close 
to unity to ensure that we are still in the weakly inter- 
acting regime. For very small momenta (we analyzed 
momenta down to k = 0.045/c G for k c /k G 45), the 
real part of the normal Matsubara Green's function in- 
dicates a resonance at u> — ck; however for such small k 
the analytic continuation no longer converged. We never 
found any trace of the Ginzburg scale k G in the proper- 
ties of the normal Green's function, which, at least in the 
weakly interacting regime, is perhaps expected since in 
the asymptotic regime the non-analytic terms are known 
to cancel out in the normal Green's functions, see the 
discussion in Appendix ID] 

In order to check convergence of the Pade approach 
we calculated Matsubara Green's functions for up to 450 
non-equidistant frequencies in the range < uj < \QQck 
for each k and investigated the stability of the results on 
increasing the number of Pade polynomials. For exam- 
ple, for k G ~ 0.2fc c and k ^ 0.1fc c we obtain excellent 
convergence and a spectral density which obeys the cor- 
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-0.6,1 



2 

co/ck 



300 



o 
< 



200 



100 




1.0 



1.1 



1.2 



co/ck 



0.08 

LLJ 

^0.04 



°6- 



0.4 

k/k„ 



0.6 



0.8 



FIG. 6: Real and imaginary parts of the normal Matsubara 
Green's function as a function of cj, calculated within the non- 
self-consistent FRG approach. The parameters are Uo = 0.8 
and /2 = 0.08, which yields k c /k G ~ 45. The Green's func- 
tion is calculated for k = 0.01fc c « 0.45fc G . The dashed verti- 
cal line shows the position of the minimum of the imaginary 
part, which coincides with the predicted resonance at u — ck. 



rect normalization 



2^ 



A(k,cu) = 1 



(4.2) 



FIG. 7: FRG results for the single-particle spectral density 
A(k,uj) as a function of the real frequency uj for different val- 
ues of k (top figure) and quasiparticle damping (bottom 
figure). The shown results are for /io = 0.15 and uo = 15. In 
this case A?g/A; c ~ 0.2 which is indicated by the labels in the 
figures. The inset in the top figure shows the quasiparticle dis- 
persion Ek which deviates at large k from linearity but is well 
described by a Bogoliubov type expression Ek = c2 ^ 2 
with renormalized velocity c = (2mVZ)~ 1 ^ 2 (black dots). 
The peaks of the spectral function correspond to (from left 
to right) k/k c = 0.2, 0.3, 0.4,0.5 and 0.6, where k c = 2mc. 
In the lower plot, black dots are extracted from the spectral 
density, while the solid line fits them as ■y^ « 0.194fc 3 /2mfc c . 



to high accuracy. The convergence properties of the Pade 
approximation become continuously worse on lowering k 
and we could not get a satisfactory analytic continuation 
in the non-perturbative regime fc <C kc- However, this 
may be simply due to the narrowing of the peak width 
for small momenta. For all momenta where we do get a 
stable analytical continuation we find a damping which 
agrees qualitatively with the expected Beliaev damping, 
even for larger interaction strengths. 

A typical shape of the spectral density function ob- 
tained from the Pade approximation is shown in Fig. [7] 
for positive u>. One clearly observes a finite peak broaden- 
ing which can be ascribed to Beliaev damping^ The ex- 
tracted spectrum of elementary excitations is always well 
fitted by a Bogoliubov-like expression Ek = \J t\ + c 2 fc 2 
and the damping of quasiparticles always reveals a fc 3 
behavior for fc — ¥ 0, in accordance with the predictions 
of the perturbative analysis, compare Eq. (|2.19[) . How- 
ever, the prefactor ao introduced in Eq. (|2.19l) should be 



replaced by a function a(pio,uo) of the relevant dimen- 
sionlcss parameters /2q and uq of the model (|2.1[) . 



(£>= 

7fe 



:2) 



a(p, ,u )k 3 



(4.3) 



For small values of both dimensionless parameters fio 
and Mo our results for the damping are very close to 
the perturbative result Eq. (|2 . 18|) . For instance, for 
flo = 0.008 and uo = 0.8 we obtain c/cq w 1.0054 and 
a(flo , uq) / c^o ~ 0.967. As both parameters increase (at 
fixed UV-cutoff Ao), we find a stronger renormalization 
of both the velocity c of the Goldstone mode and of 
a(flo,uo). For instance, for jEto = 0.4 and uq — 4 we 
obtain c/co ~ 1.01 and a(jlo,Uo)/ao ~ 0.915, while the 
results for £lq = 0.15 and uo — 15 are c/co ~ 0.669 and 
a(flo,uo)/aQ w 0.526. The first examples correspond to 
the weakly interacting regime, where many-body renor- 
malizations are expected not to play an important role, 
while the last choice of the initial conditions corresponds 
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to a more strongly interacting Bose gas where deviations 
from standard perturbation theory are already notice- 
able. 



V. SUMMARY AND CONCLUSION 

In conclusion, we have shown how the FRG formalism 
can be applied to calculate the single-particle spectral 
density of interacting bosons in dimensions D < 3. We 
here concentrated on the zero-temperature limit, but it is 
straightforward to extend the analysis also to finite tem- 
peratures and even to the critical temperature^ The 
FRG seems to be the only presently available analytical 
technique which is capable of bridging the gap from the 
non-perturbative regime at very small frequencies and 
momenta, where the self-energies become non-analytic, 
and the intermediate to large momentum regime which 
is well described by the standard Beliaev theory>2. In addi- 
tion to satisfying all symmetry constraints, our approach 
respects the Nepomnyashchy relation for the anomalous 
self-energy^ and yields the correct asymptotic structure 
of the Green's function as first derived by Gavoret and 
Nozieres4£ Moreover, our approach not only reproduces 
the non-analytic behavior of the self-energies for scales 
smaller than the generalized Ginzburg scale kc, but it 
also yields the standard results for the T-matrix renor- 
malization in D = 2 and D — 3 at momentum scales 
larger than k c — 2mc, see Appendix B. Our approach 
therefore provides a single unified framework to describe 
the crossover from the long-wavelength Goldstone mode 
regime to the quasi-free boson regime where the energy 
dispersion is quadratic. The correct momentum depen- 
dence of the Beliaev damping, including its prefactor, can 
be extracted from the numerical analysis of the our trun- 
cated FRG flow equations. While we presented numeri- 
cal results for the single-particle spectral density only for 
D = 2, it is straightforward to apply the technique also 
in three dimensions. 

The explicit results for the spectral function pre- 
sented in this work were obtained within a non-self- 
consistent solution of the truncated FRG flow equations 
(I3.26al3.26b[) for the self-energies. Although the deriva- 
tive expansion has served as a guide for the derivation 
of Eqs. (|3.26al3.26b|) . our truncation scheme transcends 
the derivative expansion because on the right-hand side 
of our truncated FRG flow equations all powers of mo- 
menta and frequencies are retained. See Ref. [21] for a 
similar truncation strategy in the symmetry broken phase 
of classical 4 -theory. 

Clearly, our truncation strategy will eventually fail at 
strong coupling. However, in principle the FRG ap- 
proach, which is a non-perturbative RG technique, is well 
positioned also for the strong coupling regime. In fact, 
our truncated FRG flow equations (|3.20l3.23al3.23b|) 
takes a two-body density-density interaction of arbitrary 
strength into account, but neglects all higher order (irre- 
ducible) correlations involving three or more particles. A 



solution of these flow equations, or an improved approxi- 
mative solution, is nonetheless expected to give valuable 
insights to strongly coupled bosons. 
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Appendix A: Elements of Nepomnyashchy theory 

Gavoret and Nozieres showed that both normal and 
anomalous propagators have the same low-energy asymp- 
totic behavior with merely a different sign^i 

G N (K) ~ —G A (K) ~ — ™f K^O, (Al) 

where c is the velocity of the Goldstone modes, p° is 
the condensate density and p denotes the boson density. 
Note that expression (|A1I) is an exact asymptotic result. 
However, as was mentioned by Gavoret and Nozieres 
themselves, these asymptotic formulas were initially ob- 
tained by omitting divergent terms which appear in the 
diagrams for both self-energies. Gavoret and Nozieres 
further assumed that S (0) 7^ 0, which contradicts the 
Nepomnyashchy relation.- Eq. (|A1|) was later rederived 
in a formalism which is free of IR divergences and is con- 
sistent with Yj A {Q) — 0, see Ref. @. We sketch here 
briefly the arguments which lead to the result E^(0) = 0. 

The anomalous self-energy can be represented as a 
sum of one-particle irreducible skeleton diagrams, some 
of which are regular and some IR divergent, 




FIG. 8: Diagrams giving rise to the IR-divergent contribu- 
tion (|A2j) to the anomalous self-energy. Black dots denote 
the bare vertices while the shaded triangles correspond to 
fully renormalized vertices. 



Z A (K) = I1 A (K) + X A (K), (A2) 



13 



where the first term on the right-hand side contains terms Eq. (IA2|) is IR divergent and arises from the diagrams 
with a single closed Green's function. These diagrams shown in Fig. [5] which have the analytic form 
are regular for K — > and D > 1. The second term in 



S s l g (X) = 2u Aa J {r^ 2) (Q,K, Q - K)G N (Q)G N (Q - K) + Y^{Q - K; Q, -K)G N (Q)G A (K - Q)| , 



(A3) 



r 



where and T^ - 3 ) are the exact irreducible vertices 

with two in- and one outgoing, and with three ingoing 
legs, respectively. The first factor u\ in Eq. (|A3|) is 

(1 2} 

the value of the bare vertex r Ao ' . In the limit K = and 
for small internal momenta Q we employ the asymptotic 
formulas (|A1|) and obtain 



e a (o) = [r (1 < 2) (o,o,o)-r (0 < 3) (0,0,0) 



IV 

X K D J dqq - 1 J 



du 2u Ao 1 



2tt {lu 2 + c 2 q 2 ) 2 



, (A4) 



where the upper limit of the integration over momenta q 
is irrelevant since we are only interested in the IR behav- 
ior of the integral. The expression in the brackets on the 
right-hand side of equation (IA4|) can be further simplified 
by using the Ward identity^ 



r (1 ' 2) (o,o,o)-r (0 - 3) (o,o,o) 



E A (0) 



(A5) 



which, upon insertion of Eq. (|A5[) into Eq. (1A4[) . yields 
an exact self-consistent equation for the anomalous self- 
energy. Since the integral on the right-hand side of equa- 
tion (|A4[) diverges for D < 3, the only finite solution of 
this equation is 



E A (0) = 0. 



(A6) 



Appendix B: Analytic approach to the z — 2 regime 
of the derivative expansion 

Here we want to investigate the flow of the interaction 
constant u A in the regime where z A = 2, i.e. the ini- 
tial part of the renormalization group flow. This regime 
yields results which are analogous to the renormalization 
of the T-matrix in both D = 2 and D = 3. The two- 
dimensional situation is especially interesting since the 
scattering length vanishes in D = 2. We work out below 
the relevant scale which replaces the scattering length 
and recover the results first discussed by Schick 3 ^ who 
investigated the hard-core Bose gas in two dimensions 
using diagrammatic techniques. In D — 2 the vanishing 



of u\ for A — > in the initial z = 2 regime is logarith- 
mic in A/Ao, which reflects the vanishing of the s-wave 
scattering length. 40 In the z — 2 regime, it suffices to 
investigate the flow of the interaction, while setting all 
remaining coupling parameters of the derivative expan- 
sion to their initial values. Performing the integration 
over frequencies we obtain from Eq. (|3.35b[) 



d A u A 



K D+1 K D e 3 



5^A A 



2mD 



E 5 



^,(Bi) 



where E A = \Aa( £ a + -^a) and A a = p A u A - In the 
z = 2 regime, the Bogoliubov spectrum approaches the 
free particle dispersion, i. e. e A 3> A A . In this case we 
may simplify Eq. (|F31[) as follows 



^a 1 ~ -2m{K D /D)A D - 3 . (B2) 
For D > 2 the solution of Eq. (|B2I) is of the form 



u A 



'■ A 



2mu Ao K D _ ^ D _ 2 ^ 



(B3) 



1 + 



D(D - 2) 



For A — > the interaction flows to a finite value which 
for u Ao Aq~ 2 ^> 1 is given by 



~ 2mu A K D n „ 
1 + . A « U A°~ 2 



(D-2)a 
2to(2tt) 



D-2 



(B4) 



D(D - 2) 



where we related the UV-cutoff A rj 2-k ja to the inverse 
s— wave scattering length a. In D = 3 we obtain 



3tt 
2m 



(B5) 



which is smaller than the scattering theory prediction for 
the T— matrix 8na/m (see for instance Ref. [|j]). How- 
ever, this can be accounted for by re-adjusting the UV- 
cutoff to A = 37r/8a. According to Eq. (IB3|) . the s-wave 
scattering length is finite in all dimensions above two. 
In D = 2, the solution of Eq. (|B2|) has the form 



it 



mUf. ( An 

1 + ^ In 

2tt U 



(B6) 
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In the regime 



1 > 



4 A 



2tt 



(B7) 



the interaction does not flow and the Bogoliubov theory 
remains correct. The opposite case 



1 < 



2tt 



(B8) 



corresponds to the limit ol hard-core bosons, where the 
flow does not depend on the initial value of the interac- 
tion. In this case we may write 



2tt 



rain (Aa 



(B9) 



which reproduces the result obtained by Fisher and 
Hohenberg4S If we now define the crossover scale where 
the scaling behavior of the Bogoliubov spectrum changes 
as e A ps 2p°u A , we find, setting up to logarithmic accu- 
racy ln(A /A c ) 2 fa 1, 



A 



16-7T/9. 



(BIO) 



Hence, the corresponding energy crossover scale is given 
by 



A Ac pa pu Ac ps - 



8np 



rain (p°d 2 ) 



(Bll) 




which is Schick's crossover scale for hard-core bosons^ 
where d is of the order of the diameter of the hard-core 
disc and we chose A = ^16n/d 2 . The velocity of the 
Goldstone mode in the hard-core regime becomes 



(B12) 



Thus, the FRG flow equations qualitatively describe the 
hard-core limit. 



Appendix C: Asymptotic behavior of the derivative 
expansion 

Here we derive the asymptotic A — > behavior of the 
flow equations obtained within the derivative expansion, 
Eqs. (|3.35al I3.35b[ I3.39aj|33^c|) . and derive from them 
the Ginzburg scale k& which is characteristic for this 
regime. This regime differs from the one discussed in Ap- 
pendix [B] in that now the dynamical exponent is z\ = 1. 
The non-perturbative character of the Bose gas is mani- 
fest in this regime. 

For A — > we may replace all coupling parameters 
which remain finite in this limit by their fixed point val- 
ues, i.e. Z\ ss Z, Va ~ V, p A ps p°. We may further 
approximate r] A ps since rj A plays no important role 
at T = 0. We are thus left with the flows of u, and 



Y A . Assuming |w| ^5 cA, where c denotes the renormal- 
ized velocity of the Goldstone mode defined in Eq. ()3.46|) , 
the small A behavior of the denominator in the propaga- 
tors (j3~^5aj) and (|3.28b|) is^ 

Ylu? + (e A + V> 2 )(2A A + e A + V A u?) 
~ 2A A (e A + Vlo 2 ) = 2VA a (lo 2 + c 2 A 2 ). (CI) 

For small A, the leading contribution in the numerator of 
Eq. (l3~35b| arises from the term 4A A in Cg ^ defined in 
Eq. (|3.38a|) , since all remaining terms contain additional 
powers of e A oc A 2 . The integration over frequencies can 
be carried out and the flow equation of the interaction 
Eq. (|3.35b|) simplifies to 



d A u A « A 



A A D—A 



A, 



where 



AD 



3mk c KD \P' 



mk r . 



(C2) 



(C3) 



Here k c = 2mc and p s — Zp° is the superfluid density 
which at T = coincides with the density of bosons^ 
and 



A' n 



AD 

\p 



(C4) 



One has to distinguish between the D < 3 and D = 3 
cases. In the first case, the general solution of Eq. (IC2[) 
can be found in the form 



ua 



u r A e 



A f 



mk c u Q 



cA 



D 



(C5) 



where e = 3 — D. Here, we introduced the scale A c and 
ua c = u c . The scale A c is the characteristic scale for the 
boundary of the perturbative z — 2 regime and can be 
approximated as A c « k c . At this scale, ma has already 
been renormalized by the RG flow through the pertur- 
bative z = 2 regime, loosely defined by A c < A < A . 
This in general can lead to a substantial suppression of 
u c from its bare value ua - However, for weak coupling, 
for an order of magnitude estimate it is sufficient to as- 
sume u c ps ua - In any event, the flow of the interaction 
for A A c in D < 3 neither depends on the chosen 
initial interaction u c , nor on the starting value for the 
cutoff-parameter A c , but only on the fully renormalized 
quantities entering the momentum k c and the factor A' D . 
For e > 0, i. e. for D < 3, the asymptotic expression for 
ua is in this limit given by 



"A^o 



**D 
mk r . 



A' 



(C6) 



For D < 3 the crossover from the z = 2 to the z = 1 
regime is governed by the generalized Ginzburg scale kc, 
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FIG. 9: Comparison of the full numerical solution of the flow 
equation for the interaction ma as obtained from Eqs. (|3.35al 
I3.35bl I3.39all3.39cl) (solid lines) in D = 2 with the approxi- 
mate solution Eq. (|C5|) (dashed lines) for different values of 
the bare parameters. Solutions 1 are calculated with u = 0.1, 
jUo = 0.001 which implies k c /k G « 340. Solutions 2 were 
calculated with u — 0.8 and p, = 0.08 with k c /k G w 45. 
Solutions 3 are obtained with u — 15 and p = 0.15 with 
k c /kc ~ 5. The vertical dotted lines show positions of the 
Ginzburg scale for each solution. 
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which can be read off from Eq. (|C5 

_ i 

mk c u 



-.4 



D 



(C7) 



For the weakly interacting Bose gas we have k c f=a 2mco 
and p s m p° . Then, ko may be rewritten in a form similar 
to the one which was obtained by Castellani et al& and 
Kreisel et a/.,— 



(C8) 



In Fig. [9] we compare of the full solution of the flow of 
ma, obtained from Eqs. (|3.35al f3.35b[ I3.39all3.39cj) . with 
the approximate solution Eq. (|C5[) for D = 2 and dif- 
ferent initial conditions. While solutions marked with 
numbers 1 and 2 in Fig. [9] correspond to a weakly inter- 
acting system (see the legend), i.e. the ratio c/c s» 1, 
where c represents the Goldstonc mode velocity calcu- 
lated from Eq. Q3.46P and c = \J p Q u Aa /m corresponds 
to the mean-field result for the Goldstone mode velocity, 
the solution marked with the number 3 corresponds to a 
more strongly interacting system with c/c ~ 0.669. It 
is evident that the approximate solution Eq. (|C5|) always 
yields a very satisfactory approximation for A A c . 

In D = 3, the only non-trivial solution of Eq. (|C2j) 
reads 



u A 



1 + 



-4; 



(C9) 



which yields the Ginzburg scale 



A c exp 



4' 

^3 



mk„u, 



(C10) 



where A c w k c . In the weakly interacting regime, this 
expression reduces to 



k c exp 



8tt 2 P ° 
~kT 



(Cll) 



which is also consistent with the result reported in 
Refs. @,|4l|. 

Analogously, from Eq. (|3.39al) we obtain, 



A 



(C12) 



D 



Taking in account the initial condition Y c as 1, the solu- 
tion of Eq. (|C12|) reads 



A' 



mk..u r . 






A,+ 4 





in D < 3, and 



mk c u ( A, 



1 



-4^ 



A 



(C13) 



(C14) 



in D = 3. Hence, the generalized Ginzburg scale is also 
characteristic for the flow of the coupling parameter Y A . 



Appendix D: Asymptotic behavior of propagators in 
the infrared limit 



We here extend the derivative expansion analysis of 
Sec. IIII Al to include also an expansion of u\(K) to second 
order in momentum and frequency. We thus approximate 



u A (K) 



+ a A e k + Pkoj 



(Dl) 



As mentioned in Sec. IIII Al the coupling function u A {K) 
is expected to become non-analytic for A — > and this 
can be understood already partially from the scaling di- 
mensions of a \ and /?a, 



[a a J 



2 - D - za, 
4 - D - 3z A . 



(D2a) 
(D2b) 



In the Goldstone regime we have z A = 1 which for D = 2 
yields [a a] = [/3a] = —1, indicating that both a a and 
Pa diverge as A -1 for small A. This suggests that for 
small K the leading order behavior of ua^>o{K) should 
be proportional to k or \u)\, which is consistent with a 
\J u 2 + c 2 k 2 behavior which we find from the full fre- 
quency and momentum dependent calculation. In D — 3, 
the dimensional analysis predicts [a a] = [Pa] = — 2. This 
is suggestive of a logarithmic dependence on |w| and k 
which is indeed correct. 

The flow equations for a A and /3 A can be obtained from 



d A a A = 
9aPa = 



Pa 

13, 



m { d 2 
~pj A [dk 2 ^ 



£?(fe,0) 



k=0 



^ A p A +^S A (^S^(0,a;) 



Pa 



2p< 



(D3a) 
,(D3b) 
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FIG. 10: 5 A = Aq «a and /3 A = 2mA" ~ 2 /3 A , calculated 
in D — 2 for the initial values of the dimensionless chemical 
potential £t = 2m/iA ( 7 2 = 0.4 and dimensionless interaction 
u = 2mu Ao Aj ) " 2 = 4. 



where we employ a non-self-consistent approach to calcu- 
late the flows, i.e. the flows are calculated in exactly the 
same approximation as the full frequency and momentum 
dependent calculation presented in Sec. IIVI The initial 
conditions of a A and (3 A are a Ao = and /3 Aq = 0. 

The analytic expression for the flows is rather lengthy 
and can be found in Ref. [42] . The flow of the parameters 
a A and (3 A for D = 2 is shown in Fig. [TO] The divergence 
of both parameters can be understood by keeping only 
the leading order terms of their flow equations in the limit 
A — > 0. In this case, both flow equations reduce to 



d A a A oc d A (3 A oc it A A 



D-6 



For D < 3 we have from Eq. (|C6[) u A oc A 3 D and thus 



the solution of Eqs. 
havior, 



D4p has the following small A be- 



a A oc (3 A cx A 1- 



(D5) 



which agrees with the dimensional estimate for z A = 1. 
In D = 3, the small A behavior ua is proportional to 
ln(A/A ) [see Eq. ([U9)l]. so that for small A Eq. (jD4|) 
leads to 



a A oc /3 A oc / dA' 



1 



1 



(A') 3 In (A /A') A 2 ln(A /A) 

(D6) 

where we have kept only the leading term. The upper 
integration boundary is not relevant for this estimate. 
Hence, using A ~ \JuP- + c 2 k 2 , we obtain the analytical 
behavior of the physical anomalous self-energy in the low- 
energy limit in 1 < D < 3, 

Y> A (K) oc (lo 2 + c 2 k 2 Y' 2 . 

In D = 3, the leading behavior of the anomalous self- 
energy is logarithmic, 



H A {K) oc 



1 



In 



cA n 



We are now in the position to recover the low-energy 
behavior of the Green's functions. By taking Eqs. (|3.9a[ 
I3.9bl 13.171 IDlj) into account, one finds 



1 



V K {K) 



[A A +iY A u + (Z A 1 +p° A a A )e k 



G A (K) = 



+ R A {k) + {V A + P lP A )u 2 ] 
1 



(A A + p° A a A e k + p A /3 A w 2 



V A (K) = Y^ 2 + (e k /Z A + R A (k) + V A u 2 ) 
x [(ZX 1 +2p A a A )e k + R A (k) 
+{V A + 2plP A )u 2 ]. 



(D7) 
(D8) 

(D9) 



The divergence of the parameters a A and j3 A enables 



us to expand Eqs. (|D71 ID8I) in powers of the small quan- 
tity a^ 1 under the assumption (3 A j ol A — > const. For 
momenta k < A, we obtain to leading order of the ex- 
pansion, 



G N (K) = -G A (K) 



2V[u 2 + {2jnZV)- 1 k 2 } 



n 2 
mp c 



c 2 k 2 ' 



(D10) 



(D4) where p — p° /Z, and the velocity of the Goldstone mode 



is 



1 



V2mVZ 



(Dll) 



see Eq. p. 461) . Eq. (|D10j) is exactly the Gavoret and 
Nozieres result* 3 ^ given in Eq. (|A1|) . We conclude that our 
non-self-consistent approach from Sec. IIVI both correctly 
describes the non-analytic behavior of the self-energies 
and yields the correct asymptotic behavior of the prop- 
agators. Note also that in the asymptotic result (|D10|) 
the scale ka and all of the non-analytic behavior of the 
self-energies is completely absent and it remains an open 
problem whether and how the non-perturbative scale kc 
shows up in the spectrum or damping of quasiparticles. 
The scale does however emerge quite dramatically in the 
so-called longitudinal Green's function which is defined 
and discussed e.g. in Refs. [6 Ul5lfl6| . While for usual 
Bose systems the longitudinal Green's function is not ac- 
cessible to experiments, in antiferromagnetic materials, 
subject to strong magnetic fields close to the saturation 
field strength, it is possible to directly probe the spectral 
properties of the longitudinal Green's function^ which 
is related to the longitudinal spin structure factor of the 
underlying spin model. 
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